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A Kosterlitz-Thouless phase transition in the ground state of an antiferromagnetic spin-^ Heisen¬ 
berg chain with nearest and next-nearest-neighbor interactions is re-investigated from a new per¬ 
spective: A mapping of the components of the scalar product onto a set of loops is found. One 
can classify these loops according to whether any two of them can be transformed into each other 
in a continuous way (i.e., whether they have the same winding number). A finite size scaling of 
the fidelity susceptibility and geometrical phase calculated within each class of above mentioned 
loops leads to the accurate critical coupling constant value and enables one to find that the critical 
exponent v = 2.000 ± 0.001. 
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Introduction. — At quantum phase transition (QPT) 
properties of the ground state of the quantum system 
change drastically due to quantum fluctuations which are 
most clearly pronounced at zero temperature. Although 
many approaches have been proposed to examine QPTs, 
to locate critical points and to calculate values of criti¬ 
cal exponents, one still asks an important question: Is it 
possible to explore the critical behavior of a system at 
QPT by examining the change of its ground state j'l'o) 
in a critical region, especially when there is no possibility 
to identify an order parameter nor to establish a pattern 
related to a symmetry breaking? Still there exists, in this 
field, a quest for new approaches, based on scaling and 
renormalization, which enable searching and characteriz¬ 
ing QPTs. Recently, mostly due to interplay between in¬ 
formation theory and quantum many-body physics, new 
possibilities have been found for the studying of QPTs, 
see e.g.. Refs. [1, 2]. 

The model under consideration and its essential prop¬ 
erties. — An additional answer to the question raised in 
the Introduction is presented below result of reexamina¬ 
tion of a quantum phase transition in the knowi>^“— one¬ 
dimensional spin-i antiferromagnet with nearest (Ji) 
and next-nearest (J 2 ) neighbor interactions 

H = Ji E SiSi+i -f J 2 ^ SiSi+2. (1) 

i i 

The ground state of this system depends on A = J 2 /J 1 . 
Namely, for A < Ac = 0.2411 it exhibits the characteris¬ 
tics of a ID antiferromagnet with nearest-neighbor inter¬ 
actions only (spin liquid phase), i.e., it is critical with cor¬ 
relations decaying oc log(r)/r (r stands for spin-spin sep¬ 
aration). Excitations are gapless - the finite-size triplet 
gap scales oc 1/L (also with log correction), with L being 
the chain length. For A > Ac the system displays com¬ 
pletely different ground state (dimerized phase): correla¬ 
tions decay exponentially and triplet gap remains open 
in thermodynamic limit. The transition between these 
phases is known to be of Kosterlitz-Thouless type. 

Resonating Valence Bond (RVB) basis and winding 
numbers. — Examining this model in RVB basis sheds 


additional light on its critical properties: it enables one 
to find the critical value of coupling constant and the 
critical exponent iz. Recall then, briefly, the essential 
features of RVB approach^ to quantum spin-i systems: 
Matrix elements of the Hamiltonian are calculated not in 
the Ising basis but in the (complete) non-orthogonal basis 
\ck) taken from an over-complete set of singlet coverings: 

(Cfc|S,.S,|Q) = (-l)'^(±^)(Cfe|Q), (2) 

(cfc|ci) = with Af{ck,ci) being a number of 

loops arising when the coverings {ck\ and |c/) are drawn 
simultaneously on the same lattice (“transition graph”) 
and Ms= ^ stands for the number of singlets in the sys¬ 
tem. All singlets belonging to \ck) are oriented; d denotes 
the number of disoriented ones one meets while moving 
along the loop in {ck\ci) containing i and j. Finally, -|-| 
is taken if there is even number of dimers between i to j, 
— I in the opposite case. In order to find the ground state 
of Hamiltonian m one solves the generalized eigenprob- 
lem: i?|'I'o) = i?oC|'I'o), with C being the matrix formed 
from scalar products (ck\ci). 

Although, due to nonorthogonality, this procedure is 
a bit more burdensome from calculational point of view, 
it has a great advantage because of its clear meaning, 
especially when one uses periodic boundary conditions 
and maps the periodicity of the Hamiltonian onto a cir¬ 
cle. Then one can characterize the ground state |'ko(-^))) 
also by some topological properties in the following way. 
The scalar product, (4'o(A)|'ko(A)), can be interpreted as 
a set of loop coverings (transition graphs) of the lattice. 
One can define the winding number— for a loop (in the 
plane, around a center of the circle which results from 
periodic boundary conditions) as an integer represent¬ 
ing the number of times one encircles the center while 
moving along this loop. The winding number for a given 
transition graph is the largest winding number from the 
the winding numbers assigned to all loops forming this 
transition graph. An equivalent definitiorii^ states that 
each loop from the transition graph can carry a ,,lattice 
flux”, which can be directly related to a winding number 


2 



FIG. 1: Examples of representing scalar products {ck\ci) 
of two basis functions (cfe| (red) and \ci) (blue) in a 
12-spin system. The “transition graph” on the left is 
non contractible (contains 3 loops, absolute value of its 
winding number equals to 1, topological sector Wl), 
whereas on the right is contractible (contains 6 loops, 
its winding number equals to 0, topological sector WO). 




FIG. 2: Scalar product ('ho(^)l'I'o(-^)) in 24-spin system 
split between topological sectors WO and Wl (left). 
The ground state energy Eq{\) divided into sectors WO 
and Wl in systems up to 24 spins (right). 


of this transition graph. In the case under consideration 
there exist two possibilities: absolute value of winding 
number equals to 0 (WO) or 1 (Wl), see examples in 
Fig. m Gonsequently, each transition graph {ck\ci) may 
be classified according to its winding number. Two tran¬ 
sition graphs with the same winding numbers belong to 
the same homotopy class, i.e., they can be transformed 
into each other in a continuous way. 

The scalar product, the ground state energy and the 
derivatives of the scalar product in topological sectors. — 
Let us take into account that the ground state is normal¬ 
ized: (4'o(A)|4'o(A)) = 1. Differentiating with respect to 
A one obtains (9 a4'o|^o) = —(4'o|9a4'o). Expressing the 
ground state in the RVB basis \ci) 

|4'o) = ^ai(A)|c,), (3) 

i 

and taking into account that for the system under con¬ 
sideration all coefficients ai{X) may be chosen to be 
real numbers, leads to the conclusion that (9 a'I'o|'Lo) = 
('I'o|i9a^'o) = 0- (Another possibility is that it is purely 
imaginary. We examine the first possibility, the report 
on the second one — which leads to the same conclusions 
will be presented elsewhere). All the terms entering the 
scalar product (SP) 

(^ol^'o) = (4) 

can be classified according to their contractibility (WO 
or Wl). This is also true for derivatives (9a'Lo|'Lo) and 
(9 a'I'o|3a'I'o) since to calculate them one has to replace 
terms aiaj in Eq. (g]) by ^-cxj and respectively. 

The scalar product, ('I'o(A)|4'o(A)), calculated for 24 spin 

system, split into two topological sectors, WO and Wl is 
presented in Fig. [2j Note that contributions (from sec¬ 
tors WO and Wl) are symmetrical about the the straight 
line i and that for A = A (i.e., for the exact Majumdar- 
Gosh ground state) there are no non contractible loops in 
(4'o|4'o) — for finite L they decay oc 2“^/^. Let us stress 
that the probability of finding a “non contractible” term 
in the ground state monotonically tends to 0 as A —^. 


Considering the dependence of the GS energy Eq on A, 
one can also identify the contributions to the energy from 
the different topological sectors: Eq = (4'o|Ll|4'o)|^.j^p + 
= E^°{\)+E^^{\), see Fig. H (right). 
It follows that the GS energy of the system is associ¬ 
ated with two types of topologically different objects: 
contractible and non contractible (cf. vortex-anti vortex 
pairs and vortices at the Kosterlitz-Thouless phase tran¬ 
sition for classical 2D XY Heisenberg system) and there 
exists an energy gap (Ai?o(A) = E^^ — E^^) between 
non contractible objects and contractible ones. Growing 
A raises the energy of non contractible objects and total 
GS energy and lowers the energy of contractible objects. 

One can recall here the energy-entropy (Peierls) ar¬ 
gument pointing that the interplay between contractible 
and non contractible objects leads to a topological phase 
transition in the system under consideration. Let us esti¬ 
mate how does the the energy, AEq{X) needed to create 
a non contractible object from a contractible one (i.e., 
to pass from WO to Wl sector) depends on L. A stands 
for a pseudocritical value of A which is calculated for 
a inflection point in the dependence E^^{X) (Note that 
A —>■ Ac for L —oo). The difference, AEq{X), increases 
linearly with L: AEq{X) oc 0.3088L. The creation of a 
non contractible covering may be done approximately in 
iC'2/2 ways {Cl /2 = (l/ 2 )!(l/ 2 +i)! )■ One has 

0.3088L-Aclog(ic2/2) = 0, (5) 

and eventually, for large L, Ac = 0.223 (which should be 
compared with 0.241). 

Fidelity susceptibility and estimation of the critical 
properties of the system. — One can expect that the 
inflection points present in the scalar product (4'o|'I'o) 
dependence on A should result in extremes of its deriva¬ 
tives: (cIa'I'oI'I'o) and (4'o|i9a'I'o) in both topological sec¬ 
tors. Indeed, numerical results show that the derivative 
(calculated numerically with an accuracy 10“®) of SP 
with respect to A exhibits a peak (see Fig. [3] - top left) 
located at A*(L) > Ac. 

Having determined the dependencies of the derivatives 
of the scalar product on A, one can calculate the fidelity 
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FIG. 3: The absolute value of the scalar product 
derivative |('I'o|cI'I'o)| with respect to A in sectors WO 
and W1 (top left) and and the fidelity susceptibility per 
spin, (top right) in systems up to L=24 

spins as a function of the parameter A = J 2 /J 1 . The 
fidelity susceptibility peak heights as a function of L in 
the WO sector (bottom left) and as a function of log(L) 
in the W1 sector (bottom right). 


susceptibility (per spin), 

X/(A) = = xr w + xr'(A), (6) 

where and xY^ stand for the fidelity susceptibility 
components from sector WO and Wl, see Fig. [3]- top 
right. One can observe quite different scaling of maxi¬ 
mum of X/(A) with respect to L in both topological sec¬ 
tors: xY^ logarithmically divergent (oc 0.155log(L)), 
xY^ finite systems approaches its infinite-size value 
bo = 0.0459 ± 0.0473 with a logarithmic size correction 
oc ~ 0.380, 62 = 6.82. Thus, for L —00 

X/(A) is logarithmically divergent. This differs from 
the resultsiS obtained for ID XXZ Heisenberg spin-A 
system with open boundary conditions (OBC), pointing 
that fidelity susceptibility does not diverge at Kosterlitz- 
Thouless QPT in ID system and numerical results were 
interpreted as logarithmic finite-size scaling corrections. 
Let us stress that if we would examine the scaling of x/ 
with respect to L without its splitting into topological 
sectors, we would have received also the finite value of 
Xf for L —>■ cx) with logarithmic size correction (likeA^ 
for OBC). This correction, 6 i/log(L/ 62 ), is still present 
in the WO sector . This means that while examining 
topological phase transitions it is better to use periodic 
boundary conditions since even for large systems with 
OBC it is very easy to confuse a logarithmic correction 
with a logarithmic growth of the peak. 

In order to find Ac and the critical exponent v (for the 
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FIG. 4: Fidelity susceptibility in Wl sector before (left) 
and after rescaling (right) for systems with L =16-24 
spins. The scaling collapse leads to the optimal values 
of Ac and exponent v. The errors were estimated by 
Hnding the collapse several times taking the numerical 
data with Gaussian noise with standard deviation equal 
to the accuracy of numerical differentiation. 


correlation length) one takes into account an argument 
that x^^(A), after appropriate rescaling of arguments 
and function values, for different system sizes should col¬ 
lapse onto the same curve. Let us assumeA d^d^ that the 
fidelity susceptibility in the sector Wl, at the critical 
point, scales as 

xr(A) = d.2/'^-"‘i>(i'/"(A-Ac)), (7) 

with $ being a scaling function. The data from the left 
side of Fig. 0] has been plotted on the right side using 
rescaled values of of argument L^/^(A— Ac) and function 
Xf^^(A)L“^/‘"+^ The collapse occurs for Ac = 0.2412 ± 
0.0002 and 1 / = 2.000 ± 0 . 001 . 

Relation to geometric phase. — It was shown that 
the derivative of the geometric phase calculated for 
the ground state of the quantum ID XY model obeys 
scaling behavior in the vicinity of a quantum phase 
transitio n^^d^ . In what follows a numerical argument 
is presented that the same is true for the the system 
under consideration and enables one to find its critical 
properties. Let us number the spins consecutively along 
the chain in such a way that odd (even) numbers are as¬ 
sign always to n.n.n. Now rotate spins on odd and even 
positions around the z direction by angles (j) and — 0 , 
respectively. This defines a new, (/)-dependent basis in 
which the spectrum of the Hamiltonian (IT|) remains un¬ 
altered. Under this transformation each singlet entering 
to a given basis function transforms as follows 

Is) = ^(in) - i;t)) ^ k') = - e-*^iit)), 

( 8 ) 

and the geometric phase, J{s'\^)d(j), for two spins-A 
coupled antiferromagnetically in their ground state, due 
to equal contributions (with opposite sign) from the term 

(e* 21 'l'l) and the term e“®^|4,t)) equals to 0. The same 
is true for any finite system, but one can overcome this 
difficulty, by accumulating the geometric phase for only 
one term from Eq. m and similarly for larger systems. 
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FIG. 5: The maximum value of the derivative at 

dX 

the pseudocritical point A* (from Fig. |3]- top left), as a 
function of lattice size for the considered here system - 
squares and, for comparison, for the ID spin-i XY 
system - diamonds (left). The value of the function 
B{X,X*,L) defined by Eq. ([9|) versus L^/‘"(A — A*) for 
system sizes L =16-24. As expected from the finite size 
scaling ansatz the data for different system sizes 
collapse on a single curve for v = 2.003 (right). The 
error is estimated as explained in the caption of Fig. 2] 


Eventually, the geometric phase, /3g, per singlet, for the 
ground state |4'o(^i>)), accumulated by varying the angle 
(/) from 0 to TT is proportional to the —*7r('I'o|4'o) and 
seems to be not A-dependent one. The absolute value of 
its derivative, however, calculated in sectors WO and W1 
displays a well-marked peak at X*{L) > Ac, which shifts 
towards Ac with increasing L, see Fig. [3] (top left). The 
value of ^ calculated for A* (L) diverges logarithmically 
with increasing L: -§^l3g oc 2.272 logL, see Fig. [5] (left). 
Similar dependence was reportedly for the ground state 
of the quantum XY chain, shown also in Fig. 0 (left). The 
logarithmic divergence suggestsl^ that it is also possible 
to extract a critical exponent v from the scaling of the 


function 


7?(A,A ,T) = (l — e ^ ^ j 

(xL^/^XiL)- X*iL)), 


(9) 


X*{L) stands for the peak position for a given L (Fig. [3] 
- top left) and l3g{X{L)) - for the geometric phase in the 
vicinity of Ac for a given L . All the data for systems 
with different L collapse onto a single curve, see Fig. [5] 
(right), after taking Ac = 0.2411 for v = 2.003 ± 0.009. 

Summary. — In this paper an example is given how 
to determine the values of critical couplings and critical 
exponent in a frustrated ID spin system at topological 
quantum phase transition. The approach was based on 
finite size scaling of fidelity susceptibility and geometric 
phase; both quantities were calculated for different topo¬ 
logical sectors. Topological sectors were defined taking 
into account contractibility of loops, existence of which 
is a natural consequence of the use of the RVB basis. 
This enabled to interpret the transition as a topological 
one and to find critical couplings and critical exponent v 
for relatively small systems. We hope that the presented 
results will stimulate further exploration of critical phe¬ 
nomena in systems in which there is no possibility to 
define an order parameter. 
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